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The elastic properties of hep ''He samples have been shown to display various anomalies. As 
the temperature is lowered below ~ 0.2 K, the elastic shear modulus appears to stiffen and the 
moment of inertia appears to drop in a concomitant manner. The former phenomenon is taken as 
evidence for the pinning of dislocations, the latter, for the appearance of supersolidity. The close 
relationship between these two observables is studied within the framework of classical deformable- 
body mechanics. A model based on the formation by plastic flow of extremely soft quasi-planar, 
inter-connected layers of dislocations is solved analytically and numerically. This model relates 
quantitatively the change in moment of inertia to the drop in elastic constant and can account for 
most experimental observations. Other situations, in which its relevance may seem more doubtful, 
are discussed. 

PACS numbers: 67.80-s, 61.72.Hh, 62.40.+i, 67.90.+Z 
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I. INTRODUCTION 

Helium makes an intriguing solid. Both '^He and ■^He 
isotopes crystallise from the liquid at absolute zero un- 
der a pressure of 24.5 and 34.5 bars respectively. They 
form hexagonal close-packed (hep) crystals with similar 
elastic properties. These crystals are very soft, owing to 
their small densities and weak atomic interactions. The 
paucity of slip directions causes crystals with the hep 
structure to be quite brittle: they fracture very readily. 
This combination of softness and brittleness in the helium 
crystals make them very prone to plastic deformation and 
the formation of dislocation lines. 

Because of the fast boson exchange in solid ''He and the 
presence of defects, the possibility that a Bose-Einstein 
condensate would form within the crystalline lattice be- 
low a certain temperature was raised by a number of au- 
thors, starting with Penrose and Onsager in the fifties.' 
These authors argued that superfluid coherence, or off- 
diagonal long-range order (ODLRO), would not occur in 
an ideally perfect crystal but possibly could in a distorted 
lattice. Although the proof that they gave was criticised 
by others, '^' it marked the beginning of a long-lasting 
search, both theoretical and experimental, for features 
that could reveal the existence of such a "supersolid" state 
in a suitably disordered crystal. 

This search received strong renewed impetus after the 
observation by Kim and Chan (see the reviews [')]) of an 
anomaly in the rotational inertia of ^He solid samples as 
seen as a period shift in high-quality-factor torsional os- 
cillators (TO). The increase in the period of the oscillator 
resonance below a temperature of ^ 0.2 K, now observed 
by many groups, is taken to signal the decoupling of part 
of the helium mass from the motion of the oscillator bob. 
This effect, first discussed by Leggett ' and called in the 
present context non-classical rotational inertia (NCRI), 
occurs in a number of TO experiments with widely dif- 
ferent sizes and geometries including very confined ones 
such as those of Vycor or sintered gold. 



In the framework of the time-honoured two-fluid model 
for superfluidity, " such an observation would signal the 
appearance of a superfluid-like fraction in the solid. Such 
a "condensate" would settle to rest and decouple from the 
oscillator walls as the temperature is lowered, reducing 
the moment of inertia. This interpretation is born out 
by the fact that, if the oscillator geometry is modifled 
by a partition blocking the closed loop along which the 
superflow is supposed to take place, '■''^ the effect disap- 
pears. Also, NCRI is not observed when the oscillator 
is flUed with ^He instead of '^He, ''" which constitutes a 
strong hint that quantum statistics plays a fundamental 
role. 

The TO measurements do seem to suggest that some 
form of superfluid behaviour occurs in solid ^He below 
0.1^0.2 K but other unambiguous manifestations of the 
existence of a true superfluid component, such as a non- 
dissipative DC-flow,'"'"''* or a persistent current, a second 
sound or fourth sound'''' mode, the fountain effect, the 
signature of a BEC condensate on neutron'^''"'''' or X-ray 
diffractograms '" are still lacking in spite of the efforts 
and ingenuity of many research groups. 

Shear modulus measurements in solid helium at low 
temperature provide another class of anomalous elastic 
properties. These measurements span many years, start- 
ing with the early work of Wanner et al. -" soon followed 
by others.-'" "' They have recently been extended to the 
same range of temperatures and "^He impurity concen- 
trations as the TO experiments. " A marked decrease 
in the shear modulus G takes place in most samples of 
hep '''He upon warming from T ~ 0. The magnitude of 
the softening varies from sample to sample, depending in 
particular on the ^He impurity content x^ and the cool- 
ing history. The drop in G can be spectacular, down to 
less than 20 % from the T ~ value.-'i.25,2N,29 

Day and Beamish^'' have argued that the T and 0:3 de- 
pendence of G were mimicking closely those of the period 
shifts in TO experiments. In fact, the striking similar- 
ities between the shear modulus and the TO resonance 



frequency drops make it hard not to believe that the two 
phenomena are somehow related. Direct experimental 
studies of this possible connection have led to diverging 
conclusions. ''""'^'^ 

This article'^^ outlines one possible such link between 
those two different mechanical properties of hep ''He. It 
differs from similar attempts by other workers ^ be- 
cause it recognises from the start that the large drop in G 
requires the bunching of dislocation lines into extended 
quasi-planar highly deformable sheets, as described in 
Sec. (II). The consequences of these assumed defects are 
derived analytically for the shear modulus drop and for 
the apparent change of inertia in Sec. (III). These two 
quantities can thus be directly linked to one another. Nu- 
merical values are derived in Sec. (IV), where it is shown 
that this simple model may account quantitatively for a 
number of experimental observations. The model does 
not explain readily certain classes of experiments, no- 
tably those in confined geometries, and also the absence 
of inertia anomaly in solid '^He. Some speculations on 
these topics are offered in Sec. (V). 



II. SOFT LAYER MODEL 
II. 1. Planar layers of dislocations 

As mentioned above, the two helium crystals are very 
soft: the longitudinal and transverse sound velocities, 
in the 200 to 500 m/s range, are low. They are also 
very fragile: the yield strength is of the order of 0.2~0.5 
bar,'^'''' or very much less depending on the experi- 
mental conditions.'" Dislocations appear readily under 
very weak mechanical perturbations or thermally induced 
stress. Plastic flow takes place during the formation and 
subsequent cooling of the solid helium sample. In the 
process, dislocations form and migrate. 

Early sound propagation measurements in hep '^He in 
the 5-50 MHz frequency range '.-'.^i.^- have revealed an 
anomalous temperature dependence below 1 K of the 
longitudinal sound velocity. This anomaly has been at- 
tributed to the unpinning of the dislocation lines as the 
trapped '^He impurities escape from the dislocation cores 
by thermal activation. This interpretation is well doc- 



umented through the work of many authors" 



and. 



more recently, by Syshchenko et al. ' ' . 

The analysis of the high-frequency sound propagation 
measurements yields typical values ^ 10^ cm^'^ for the 
density of dislocation lines A and 5 • lO^"' cm for the 
average distance between the nodes of the dislocation 
network Ln, assumed random and homogeneously dis- 
tributed. The dimensionless quantity Ai^ is thus found 
of the order of ^ 0.25. 

Lower frequency measurements"^'""'" ' interpreted in 
the same manner with the help of the following relation 
for the change of the effective shear elastic constant Goff 



relative to the T — Q value G ,"'"' 

G/Gcff = 1 + 24(1 - u)nKLlln'^ , (1) 

give, assuming a value of 0.3 for the Poisson ratio i^ and 
with the highest value for the orientation parameter il < 
1/2, much larger values of the quantity AL^. The shear 
modulus measurements by M.A. Paalanen, D.J. Bishop, 
and H.W. Dail ' ' were carried out at a low frequency of 
331 Hz and lead to a value of AL^ > 1.0 to 2.5 depending 
on samples. More recent measurements"*^'^'*'"''''''"' have 
confirmed these results. A softening of 86 % has been 
observed in an ultra-pure monocrystalline ^He sample by 
Rojas et al. at frequencies in the 10-20 kHz range. In 
this extreme situation, the quantity AL^ would exceed 
20 using the same values for v and O as above. 

These values of AL^, obtained at long wavelengths, are 
much larger than the upper limit for a dense hexagonal 
network of dislocations, which espouses the underlying 
lattice symmetry. As shown in the Appendix, this limit 
is l/v2 for an ideal hep network. The corresponding up- 
per limit of G/GeS as given by Eq.(l) is 1.2, which falls 
short of observations^'^'^"''-''"'-'' by a wide margin: edge 
dislocation lines escape from their preferential homoge- 
neous hexagonal network structure and become quite ex- 
tended. 

This anomaly clearly points towards the formation of 
inhomogeneous dislocation structures. The dislocation 
lines collect in dense arrays, such as the mosaic struc- 
ture that form along the boundaries between grains with 
slightly misaligned lattice vectors,"'' or, more generally, 
in extended planar structures. This rearrangement takes 
place during the formation of the hep ''He samples and 
under thermal stress during cool-down. 

It has been shown by numerical simulations of dislo- 
cation dynamics, notably by Amadeo and Ghoniem,"'^'""* 
that dislocations collect into different planar structures 
according to different applied perturbations. Planar ar- 
rays composed of sets of dislocation dipoles lying in 
planes containing the direction of the critical resolved 
shear stress form under monotonic stress conditions. 
Other types of structures, slip bands of parallel dislo- 
cation lines or dislocation cells, may appear under cyclic 
perturbations, provided, e.g., by mechanical vibrations. 
These planar defects have been observed in a number 
of metallurgical samples. Their phenomenology is well 
documented, as reviewed, e.g., by Takeuchi and Argon '" 
and others. Such dislocation substructures have also 
been observed in hep ^He by X-ray topography by Iwasa 
et al. '^ and by transmission Laue diffraction by Bossy 
et al. '" . 

These defect structures are thicker than the Franck 
networks that separate two grain boundaries of low-tilt 
angle. They are quite different from the random net- 
work assumed in the Granato-Liicke model,"'"' as already 
mentioned. They can be viewed as resulting from the 
propagation of dislocation pileups under thermal stress 
in a way similar to the formation of cracks in usual hep 
metals. ' ' Solid helium exists only under positive pressure 



and does not actually crack. Other types of extended de- 
fects appear and enable the crystal to yield in the defor- 
mation directions imposed by the rigid wall boundaries. 
The dislocation arrays formed in such a manner are 
densely packed and have a high density of long disloca- 
tion lines; dislocations interact strongly and are organ- 
ised in extended structures of parallel lines. They become 
extremely mobile when unadorned of the '^He impurities 
that pin them to the lattice at T < 0.2 K and when 
unhampered by thermally excited phonons that prevail 
at r > 0.8 K. Following the same line of reasoning that 
leads to Eq.(l), the resulting large values of AL^ lead 
to very soft and easily deformable layers. These layers 
separate regions with depleted dislocation densities but 
of enhanced crystalline quality in which deformation also 
occurs quite readily, at least in the directions of easy 
glide, but with the geometrical limit AL^ < 2^^/^. 



II. 2. Strain standing ^vaves: homogeneous case 

The simple model to be studied below assumes the 
existence of quasi-planar dislocation structures that fa- 
cilitate both plastic and elastic deformations. To make 
the problem easily tractable analytically, a fully-planar 
geometry is assumed: the helium sample is taken to be 
confined between two parallel plane walls located at z = 
and z = R and extending to infinity along the x and 
y axes. The deformation u induced in the sample de- 
pends on z and t only (see Fig. 1); the problem is one- 
dimensional and easily solvable. 

Shear stresses and strains are produced in the sample 
either by moving one plane, e.g., that at z = (which 
would be the transmitter in the shear modulus exper- 
iment) with respect to the z = R plane, held steady 
(which would be the receiver). Torsional oscillator ex- 
periments are mimicked by moving both bounding walls 
in unison, letting the sample inertia develop internal 
stresses. 

If the helium sample is homogeneous with a density 
p and a shear modulus G independent of position and 
time (no visco-elastic effect, no internal structure), the 
deformation u(z, t) obeys the following partial differential 
equation: 



a u „o u 



(2) 



This equation describes the propagation of transverse 
plane waves with dispersion relation lu'^ = c^k^ and 
c^ = G/p. The harmonic solution of Eq.(2) at frequency 
uj/2t: is the sum of two counter-propagating waves: 

m(z, t) = (uo+e^**--^ + wo-e'*^^) e^"* = u{z)e"^* . (3) 

The constants of integration for shear measurements Ug_^ 
and Ug_ are then given by 



(4) 



and the solution of Eq. (2) for the deformation as a func- 
tion of z can then be expressed under the following form: 



u(z) 



sinfc(i? — z) 
' sin kR 



(5) 



The stress in the solid is derived from the deformation, 
still disregarding the time dependence: 

du cosk{R- z) 

<^{z) = G — == -ku„G ^— - — . (6) 

az sin kR 

The force per unit area acting on the receiver is the op- 
posite of that acting on the body, namely the internal 

stress: 



Ft, 



MR) 



kug G 
sin kR 



(7) 



so that the measured effective shear elastic modulus Gog 
is such that: 



G 

Gctf 



u^G_ 
R Fn 



sin kR 
kR 



2p2 



1- 



poj^'R 
6G 



(8) 



Equation (8) describes the change of the effective shear 
modulus at finite frequency due to the elastic response 
of the body. In the limit uj — > 0, GcS reduces to G. For 
kR — n, the body is set into resonance and, as damping 
has been neglected, the effective shear modulus diverges. 
Higher frequency modes are not considered here. 

For torsional oscillator measurements this elastic re- 
sponse of the body becomes the dominant effect. In these 
experiments the two walls at z = and z — R are set 
into identical motion u^e*'^*. The solution to Eq.(2) that 
satisfies the boundary conditions 

w(z)|^=o ='«(z)U=fl, = Mo (9) 

can be written with the help of the following relations 



(10) 



■2,—ikR 



'.liikR 



2i—ikR 



1 - e"^'-'< 1 
In particular, the stress a{z) is found to be: 

cos k{R — z) — cos kz 



az 



-kGun 



sinkR 



(11) 



The quantity actually measured in the TO type of exper- 
iments is the back-action of the sample on the measuring 
device, namely the total force Fx + ^r exerted by the 
solid helium on both walls. This force is expressed, per 
unit area, by 



Fx -h Fr = cr(0) - a{R) = 2fcGu, 



1 — cos kR 
sin kR 



pRijj u„ 



1 



puj^ 
12G 



(12) 
(13) 



c,ikR 



-ikR 



c,ikR 



ikR 



The meaning of Eqs.(12) and (13) is made clear by the 
prefactor of the right-hand side Eq.(13): uf^u^ is the ac- 
celeration amplitude, pR the "bare" inertial mass Mj per 
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FIG. 1. Forms of the stationary wave in TO's experiments 
pictured as the relative displacement u(z)/uo for z varying 
from to i? for various dimensionless compliances k of the soft 
layers, taken to be equal, for R= \ cm at a frequency of one 
kHz {kR = 0.2354). The soft layers lay at zi = 0.3ii and Z2 = 
0.9-R, at the interfaces between slab 0, left, slab 1, middle, and 
slab 2, right. The dash-dot lines mark the discontinuities of 
m(«). The dash-dash line represents the elastic behaviour with 
no soft layers. 



layer Gg, she ar plane waves propagate with wave vector 
fcs = yj p/Gs- The propagation of the deformation-stress 
vector, [w(z), cr(z)], in the layer of thickness d is described 
by the transfer matrix 



M{d) 



cosfcgd (l/ZcgGg) sinfcgd 

—ksGsSmksd cos kgd 



Although the problem of finding how plane waves prop- 
agate through the stack of slabs depicted in Fig. 1 is for- 
mally solved by multiplying transfer matrices such as 
M((i), it saves a number of algebraic steps to let the 
thickness d and the modulus Gg go to zero in such a way 
that d/Gs remains constant and equal to a. The effect of 
the soft layer is then lumped into a discontinuous jump 
in the deformation proportional to the local stress, de- 
scribed, for the layer at zi, by 



u{zi + d) = u{zi) + aio-(zi) , 
a{zi + d) = ~k'^dGsu{zi) + (t{zi) , 



(14) 



In the limit d — >■ 0, the last equality expresses the con- 
tinuity of stress across the infinitely thin layer while the 
displacement experiences a discontinuity. These bound- 
ary conditions, which could have been anticipated, also 
apply to the soft layer at Z2 with slip parameter a2- In 
the following, the soft layers will be described by their 
compliances ai/R, which are such that the parameters 
Ki — aiG/R are dimensionless quantities. 



unit area and, in the square bracket, the elastic correc- 
tion at finite frequency. This "effective mass" correction 
increases with frequency up to the resonance at kR = -k 
where it becomes very large. 

The shear modulus and effective mass corrections are 
related through Eqs.(8) and (12). Taking, e.g., R = 1 cm, 
a frequency lu/27t = 1 kHz, at a density p — 0.194 g/cm^ 
for which ct = 267 m/s, the effect of shear elasticity 
on the effective mass amounts to 4.6 10^'^, which is not 
insignificant. To a drop by 20% in G corresponds an 
apparent change in the mass by ^ 10^"^. 



II.3. Soft layers 

To account for the effect of highly deformable disloca- 
tion structures, the model is extended by introducing two 
soft layers at zi and Z2 parallel to the rigid cell walls, as 
depicted in Fig. 1. The slabs of solid helium that these 
soft layers delimit have the elastic properties of the ho- 
mogeneous crystal discussed in subsection II. 2 above. 

Strain and stress are continuous functions at the in- 
terface between the slab of dislocation-free crystal and 
the soft layer. Denoting the shear modulus in the soft 



III. MODEL ANALYSIS 
III.l. Wave propagation through the sample 

With the boundary conditions, Eqs.(14), describing 
the soft layers, the propagation of the propagating and 
counter-propagating waves through the three slabs of ho- 
mogeneous crystal with shear modulus G and obeying 
no-slip boundary conditions at the walls can be found by 
straightforward algebra. 

Wave propagation in slab 0, between z = and zi as 
shown in Fig. 1, is represented by Eq.(3), which involves 
two integration constants u^^ and Ug_ , the amplitudes 
of the counter-propagating plane waves with pulsation uj 
and wavevectors ±fc. Similar solutions obtain in slab 1 
between zi and Z2, and in slab 2 between Z2 and z — R, 
involving constants u^^, Ui_, and U2+, U2- respectively. 

The amplitudes of the propagating and counter- 
propagating waves in slab 2 are linearly related to those 
in slab 0: 






(15) 



The deformation discontinuity at the soft layer at zi yields the following relations: 



J— ifczi 



Jkzi 



ui(zi) = Ui+e ""^1 +ui_e"^^i ==?i„^(l-iKifci?)e '''^^^ + u„_{l + iKikR)e 



-ikzi 



ikzi 



a{zi) 



'ikG 



("04 



-ikzi 



Jkzi 



'ikG 



("n 



'ikzi 



Ui_e 



ikzi 



Similar relations hold between 



''i+i i^i-i 



and 



Eliminating Ui+ and Mi_ leads to the following expressions for 



the coefficients Sij of the matrix that describes wave propagation through the stack of slabs 0, 1, 2: 



S^^ = (1 _ ,^kR) (1 - ^^kR) + ^k'R^ c2*fc(^2 - ^i) = s[, + z6'l, , 
5,2 = ^^kR (1 - i^kR) e2*fc^i + ^^fci? (l + ^^kR) e^'^^^ = 6[^ + %6'{, 
^21 = -i — kR { 

r 



2 1 



l + iyA;i?)c 



-2zfczi _ ^^^^ 1^^ _ ^^^^^ g-2zfcz2 = ^*^ ^ 



2 / V 2 / 4 



5*1 



(16a) 
(16b) 
(16c) 
(16d) 



The matrix A — ||(5ij||, Eq.(16), describing wave propa- 
gation in a conservative time-reversal invariant system, 
is unitary and has determinant unity: 



^\\^22 — ^\2^2\ = 1 



III. 2. Shear modulus 



(17) 



For shear modulus measurements, the no-slip condition 
at the walls reads: 



The effective shear modulus Gcff = —a{R) R/u^ follows 
readily from Eq.(19): 



(20a) 



G _ V 
G^^kR 

k^R^ 
^ ^ 6 2 

-t- Kik^Zi{R — Zi) + K2k^Z2{R — Z2) 

- KiK2k'^R {Z2~ Zi)~\ . (20b) 



'*l+«2^.2^2 



-ikR 



{Sue 



-ikR, 



+ (52ie 



U2„e 

ikR\ 



ikR. 



{5i2e 



ikR 



Sioe 



ikR\ 



relations from which the integration constants u^^ and 
Ur, can be derived: 






((5i2-<5ii)e-'fc«~((52i-(522)e 



S^^e^kR 



ikR 



(18) 



The quantities u^^ and Ug_ expressed by Eqs.(18) now 
include the effect of the soft layers and should not be 
confused with those given by Eq.(4), which do not. The 
shear stress at the receiver a{R) is given by 



a{R) ^ G 



du 

dz 



-ikGiu^.e-^^^ 



u._e''=«) 



Expressing U2+ and U2- in terms of u^^ and Ug_ using 
Eqs.(16) and (17), <j{R) takes the following simple form: 



<j(R) 



-2ikGu„ 



with 2?= (l/2i) [{5i2~5ii)' 



-ikR 



-kGujV , 

((521 - '^22)6 



(19) 



\kR-\ 



sinkR + KikRcosk{R — zi) cos kzi 

K2kRco8 k{R — Z2) coskz2 

KiK2k R sin k{z2 — zi) cos k{R — Z2) cos kzi 



For typical numerical values, such as those used for the 
graphs in Fig. 1, the first two terms of the expansion of 
l/Gcff with respect to kR, Eq.(20b), fall within 1% of 
the exact value given by Eq. (20a) . The zeroth order term 
could have been written from scratch. The correction to 
the inertial mass turns out to be less transparent and is 
considered in the next subsection. 



III. 3. Effective mass 

The derivation of the apparent inertia of the sample 
follows that given in Sec. II. 2 for the homogeneous sam- 
ple, starting from the same boundary conditions, Eqs.(9). 
Equation (10) for the integration constants is modified as: 



Si2e-^^^ + (5226''=^ - 1 



S^^Q-ikR _|_ (52ie*'=-''^ 



(21) 



the quantity u^ being the same as in Eq.(18) for the shear 
modulus case. 

Using these integration constants, the stress at each 
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FIG. 2. Inverse effective shear modulus Ges, normalised to 
the shear modulus with no soft layers G in terms of the di- 
mensionless compliance k of the soft layers, for the same pa- 
rameter values as in Fig. 1 for various values of kR. The 
dash-dash curve represents the case with R = 1 cm at a fre- 
quency of one kHz, kR = 0.2354. The symbols mark the cases 
discussed in the text, (■) for G/G^s = 9 ,k = 4.411, (•) for 
G/GcB ~ 1.6 ,K = 0.3075. The dash-dot-dash lines corre- 
spond to a reentrant branch that, for a given value of G/Ges, 
cannot be reached by adiabatic turn-on of the k^'s. 



wall takes the following form 



ct(0) = G 



a{R) = G 



dz 

ikGu^ 

du 



ikG{u^l-u^_) 



(22) 



((5i2 + 5ii)e 



-ikG 



-ikR 
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[522 + hiY 



ikR 



z=R 



^^kR^ 



-ikGu„ 



2(511^22 - ^12^21) 



{Sii - (5i2)e' 



-ikR 



{S22 - 52i)e' 



ikR 



(23) 



The total force per unit area exerted by the helium sam- 
ple on both walls is now given, instead of Eq.(12) for the 
homogeneous case, by 



Fx+Fr = a{0)-a{R) = 2ikGu„ 



2-Siie 



-ikR c AkR 
— 0226 



using again the property that \\5ij\\ has determinant 
unity. Expliciting the quantities within square brackets 
making use of Eqs.(16) and Eq.(20a), the total force on 
the walls takes the final form 

FT = Fx + FR = fcGuo^ = ^GeffAA, (24) 

U Jx 

with AA = 2 - Sne-'''^ - 522e''''^ = 2(1 - cos kR) 

-}- {ki + K2)kRsmkR~ KiK2k R sinA:(z2 — -^i) 
X I sin A:(z2 — zi) + sinfc(i? — Z2 + zi)| . 



III. 4. Stationary waveforms 

The displacement u{z) in the sample can easily be eval- 
uated using, e.g. in the inertia measurement case, the 
solution to the wave equation expressed by Eqs.(21) for 
with the following result: 



Mq"^ and u^^ 



for slab 0: u{zY^> = uo[cosfcz-|- (_B/^) sin fcz)] 



for slab 1: 
sinA:(2zi - 
cosfc(2zi — 

for slab 2: m(z)(2) 
KikR/2 cos k{z — 



u{z)^^' — Mo[cosfcz — 
{B/A)[smkz + 






KikR/2{smkz + 
KikR/2{coskz + 



Z2) coskz2 
Z2) cosfczi] 



= (uoM)[sinfc(i? - z) + sinfci? 
zi)coskzi + K2kR/2 cos k{z 
zi)cosk{z 



K.1K2K G sin k{z2 



In these expressions, 

A = {5[^- (5^2) sin kR + (J'/j - 5'1^) cos kR ^ V , 
B = l- ((5'ii + 5^2) cos kR - {S'i^ + 6'i{) sin kR . 

These waveforms can readily be evaluated numerically. 
As an example, the relative displacement u{z)/uo for 
three values of the dimensionless compliance k of the 
soft layers, taken to be equal, is shown in Fig. 1 for 
kR = 0.2354, zi = 0.3i? and Z2 = 0.9R. The discon- 
tinuities at zi and Z2 caused by these soft layers increase 
in size with the compliance, up to the point where 2? 
becomes zero and the deformation diverges. 

The next step, carried out in the following Section, 
consists in extracting the parameters Ki of the soft layers 
from the measured value of Goff and in evaluating the 
corresponding apparent change in inertia. 

IV. NUMERICAL RESULTS 

As the temperature is raised from absolute zero, k 
varies from its low temperature value, assumed to be neg- 
ligibly small because the dislocations are immobilised by 
the isotopic impurities, to its high T value. The corre- 
sponding change of G/G^tf to lowest order in the small 
parameter kR in Eq.(20b) reduces to a static correction 
to the elastic constant. The lowest order correction to the 
effective mass results from a dynamical effect of a mag- 
nitude comparable to that of the plain elastic response, 
which should be subtracted out. This difference follows 
fromEqs.(12) and (24): 



Ft 



pRuj 



'■uof 



Ki , K2— kR 

,2„, k-'R 



f\f o 1 — COS kR 

V ^ sin kR. 



I (25a) 



2 n2 



P^^O 1 + M+K2 ^ 
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zi{R-zi) 



Z2{R~Z2) 

K2 ]^2 



(25b) 




FIG. 3. Relative change in the apparent inertia, /ncri vs the 
inverse dimensionless shear modulus G/GcS, for the situation 
of Fig. (1), for various values of kR as labelled in the figure. 
The symbols (•) and (■) on the dash-dash curve for kR = 
0.2354 mark the same cases as in Fig. 2. The dash-dot-dash 
portions of the various curves correspond to reentrant regions 
that are not be reached by adiabatic turn-on of the compliance 
of the soft layers. 



Equations (25a) and (25b) show how what could 
be called the "superfluid fraction", /ncri — (Ft — 
-FtUi,k2=o)/P-^^'^o depends on the compliances Ki, 
which in turn are related to the effective shear modulus. 
The quantity pRuP'Uo has already appeared in Eq.(12) 
and stands for the force due to the acceleration of the 
inertial mass Mj = pR. These quantities hold per unit 
area. 

The full expression of the exact result, Eq.(25a), is 
fairly lengthy and not particularly transparent but eval- 
uates numerically quite readily. The outcome is discussed 
below. The lowest order correction to this effective mass, 
Eq.(25b), is second order in kR and linear in the k^'s, 
the term in kiK2, of order (kR)^, being discarded. This 
correction is either positive or equal to zero for the spe- 
cial case zi = Z2 = R/2, that is, for a vanishing dangling 
mass, and, by symmetry, vanishing local stress in slab 2. 

Equations (25a) and (25b), together with (20a) and 
(20b) , which express the NCRI fraction and the effective 
shear modulus in the presence of soft dislocation arrays, 
constitute the main result of this work. ' '■ '' 

The variation of the shear modulus in terms of the soft 
layer compliances, taken for simplicity to be equal to a 
common value k, is shown in Fig. 2 for various values of 
kR for the same sample geometry and parameter values 
as in Fig. 1. As the compliance k increases from zero, 
assumed to be its T = value, the effective shear modu- 
lus GcS decreases; the solid becomes softer, up to a point 



where G/GcS reaches a maximum: the interfaces between 
slab 1 and its neighbours becomes so soft that, although 
the dangling slab swings with increasing amplitude, the 
stress due to its motion ceases to increase. Beyond this 
point, a further increase in k would lower G/Gcs because 
the stress reflected back onto the external boundaries ef- 
fectively decreases while the displacement of slab 1 goes 
on increasing. 

It has been assumed above that the steady-state regime 
is reached adiabatically, which implies: 1) that only the 
ascending branch of G/GcS in Fig- 2 can be reached by 
adiabatic turn-on of the k^'s from zero; 2) that damping 
does not vanish entirely. If damping is introduced in 
the wave equation, Eq.(2), slab 1 would be coupled to 
its neighbours by friction in addition to shear elasticity 
and the results obtained above would be quantitatively 
different from those in Fig. 2 for very small values of Goff . 
In particular, the descending branch of G/Gcs{k) could 
not actually cross the a;-axis. 

The NCRI fraction is plotted directly in terms of the 
shear modulus in Fig. 3. From the measured overall 
change in G/GcS, which reaches values of 1.6 ' up to 9 
or more, the corresponding values of the compliance of 
the soft layers can be found from Eq.(20a). From these 
values, n = 0.31 for G/G^s = 1-6, k = 4.4 for G/GcS = 9, 
the NCRI fractions given by Eq.(25a) are 0.22 % and 0.58 
% respectively. These values depend on zi and Z2: the 
largest fncRi a-re achieved for zi '-^ 0, 2:2 ~ i?. The 
reentrant branches of the graphs in Fig. 3 correspond to 
the descending branches for large k in Fig. 2. An ac- 
curate description of these regimes where displacements 
become very large should, as already mentioned, include 
damping. They are irrelevant to the present discussion. 

Dislocations are also found in the homogeneous slabs. 
They may also induce a variation of A: = \J pjGui as G 
may also vary with temperature. The contributions to 
/ncri of the soft layers and of the dislocation network 
in the homogeneous slabs are seen in Eq.(25a) to be ad- 
ditive and their respective weights depend on how each 
contributes to G/Goff. However, the contribution of the 
network should be no more than 20%, the geometrical 
limit for hep structures, so that its effect on /ncri is less 
significant. 

The highest values of /ncri for given G/G^a and kR 
are reached for zi = 0, Z2 = ^ and varying ki while 
adjusting K2 to keep G/GcS constant. These maximum 
maxim,orum, values are plotted in Fig. 4 for kR — 0.2354 
and 0.1, and for various values of the effective shear mod- 
ulus. These quantities overtake those for homogeneous 
systems, given by Eqs.(8) and (12), up to G/GcS values 
that are much larger than those for homogeneous systems 
AL^ < 2"-'^/^, as seen in Fig. 4. Values of /ncri rang- 
ing from near zero to above one per cent can be reached 
for kR = 0.2354, i.e. in a one cm size cell at a fre- 
quency of 1 kHz for the experimentally observed values 
of G/Gaff. For a cylindrical geometry, these values are 
approximately halved. ' ' 
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eff 



FIG. 4. Maximum possible value of /ncri in terms of G/GcS 
for fc_R=0.2354 (upper curves) and 0.1 (lower curves). The 
nearly straight curves stem from the homogeneous elastic re- 
sponse, the bending curves are the pure soft layer contribu- 
tions. The values for /ncrj are extracted from graphs for 
/ncri in terms of ki as shown in the insert where they are 
marked by (a). The symbols (•) and (■) on the curve for 
kR = 0.2354 mark the same cases as in Fig. 2. 



V. DISCUSSION 

Because of the very large drop of the shear modulus 
observed at low frequencies in most samples of solid *He. 
it is surmised in this work that, instead of forming homo- 
geneous random networks, dislocations crop into organ- 
ised slip bands or quasi-planar arrays of sizable thickness. 
When the '^He impurities evaporate from the cores of the 
dislocation lines, the latter become very mobile and the 
planar arrays very soft. It is argued that a large degree 
of polycrystallinity does not suffice to obtain large G/Gctf 
and large /ncrIi ^-s shown experimentally in aerogel by 
Mulders et al. '"' . Dislocations have to organise over large 
distances in such a way that parts of the sample become 
uncoupled and acquire additional kinetic energy, thereby 
increasing the apparent inertia. 

In a hep structure, edge dislocations glide easily in 
the basal plane along three preferred crystallographic di- 
rections at 120 degrees of one another. To climb away 
from these directions, they have to change into screw or 
mixed dislocations. As shown by Suzuki and Nishioka 
this process is thermally assisted above ^0.8 K (at p = 
0.192 g/cm'^) and proceeds by quantum-tunnelling below. 
More recent and detailed theoretical considerations of the 
climb process in the quantum regime"'^'''" have led to the 
realisation that the cores of screw and edge dislocations 
could become superffuid. Climb processes would then 
become greatly enhanced, hence the concept of "super- 
climb" introduced by Kuklov and coworkers.'""'' This 



quantum-assisted climb process provides a mechanism 
in solid ''He for edge dislocations to easily move off the 
basal plane. This process lifts a constraint on dislocation 
motion. The propagation of dislocation pile-ups in the 
course of plastic deformation becomes greatly facilitated, 
as well as the formation of percolating planar defects. 

On heuristic grounds, propagation of dislocation pile- 
ups in brittle materials, such as hep helium, causes cracks 
to form and results eventually in mechanical failure. In 
solid helium, which is under positive hydrostatic pres- 
sure, cracks with voids cannot form but corresponding 
macroscopic defects with little or no crystalline order 



must appear. 



Hence the plausible appearance of con- 



nected veins imprinted by plastic ffow. These regions of 
the sample show strong spatial disorder and can possibly 
sustain off-diagonal long range order instead. '" That 
they axe found anomalously soft in a number of experi- 
ments lends credence to this possibility. 

The model based on these soft layers is easily tractable 
analytically. The calculated values for the shear modulus 
and the NCRIf fall within the range of observed values, 
barring the highest ones. This model explains read- 
ily why the NCRI and stress-strain measurements de- 
pend so strongly on the sample geometry ' and thermal 
history: '■''' even small changes in the soft layer prop- 
erties and the interconnection of the channels that they 
delimit can greatly inffuence the motion of the dangling 
masses. Homogeneous dislocation networks, besides hav- 
ing a limited effect of the shear modulus, can hardly ex- 
hibit such variability. 

Actual samples studied in the laboratory are likely to 
be more convoluted than sketched in Fig. 1. The veins 
have tortuous paths and coarse sheaths, which might 
appear to hinder motion. However, applied strains are 
small, of the order of 10~^ or less,""''"'"' and displace- 
ments are correspondingly small. The soft layers consid- 
ered here are thicker than low-tilt angle subgrain bound- 
aries, possibly in the 10 to 100 nm range. The crystal 
lattice is heavily distorted over such a thickness. The soft 
layers can be expected to be quite malleable on such a 
scale and yield easily under local stress. Being extremely 
compliant, they support plastic flow within themselves 
and accommodate departures from the ideal planar ge- 
ometry depicted in Fig. 1. The soft layers can conceivably 
also become fully fluid, '"'' or even genuinely superfluid as 
already mentioned above.'" Crystal subgrain motions on 
a sub-millimetric scale have actually been reported by 
Burns et al. in X-ray experiments using very fine coUi- 
mated beams. Similarly, the mobile features observed in 
solid ^He at higher temperatures " can be re-examined 
in the present framework; these experiments also provide 
possible clues for the existence of veins of easy deforma- 
tion. 

Maris and Balibar' ' take a quite different point of view 
to account for the observed relationship between GcS and 
/ncri- They point out that experimental TO's may lack 
sufficient structural rigidity. If the TO body deforms in 
such a way as to induce additional strain on the helium 
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sample, the NCRIf may appear larger than the intrinsic 
value. As discussed in Ref.[7i], the effect can be quite 
large. This helps in particular to understand some very 
large NCRIf values reported in the literature^' ''- that 
would not be readily explainable with the existence of 
soft layers as assumed here. From the results in Figs. 3 
and 4, /ncri remains below a few percent at the most 
for centimetre size cells, less for smaller toroidal annuli. 
But conversely, these results do not imply that, whenever 
the stiffness of the helium sample changes, an apparent 
NCRI is bound to occur; this last feature depends on 
the geometry of the soft layers and may be altogether 
absent. ' 

A number of experiments might seem to invalidate the 
present approach. The TO experiments with a blocked 
channel show a much reduced NCRIf. This is interpreted 
as the manifestation that some sort of superflow is taking 
place when flow paths are connected in a loop and not 
when the loop is broken. However, the same considera- 
tions apply to the plastic flow in connected veins, which 
also can form, or not, channels through which dangling 
masses can jiggle. 

Torsional oscillator measurements in conflned geome- 
tries, Vycor, porous gold, aerogel,. . . , do show a sizable 
NCRIf and would also appear to completely invalidate 
the present approach. If the model is applied to a single 
pore, for which kR is very small, then, indeed, the result- 
ing effect that decreases as (kR)^ will be extremely small. 
For sizes comparable to that of the soft layers, the soft 
layer model is not expected to apply, neither for shear nor 
for inertia. Pores do not appear to be filled with homo- 
geneous hep solid but with either a combination of layers 
of liquid and of bcc solid on top of 1^2 layers of amor- 
phous solid, or, for finer pores (47 A in MCM-41 and 34 
A in gelsil) , ' ' with amorphous solid only and inclusions 
of bcc-like nodules. What was assumed for softer lay- 
ers carries over to the fine pores, which present a multi- 
connected geometry with complex plastic fiow patterns. 
The conditions of existence of connected veins assumed 
at a macroscopic level are clearly fulfilled at the meso- 
scopic level of the pores so that helium, either liquid or 
amorphous, would contribute to NCRI. These questions 
deserve further consideration.^'"'' ' 

Hexagonal solid '^He is also soft but appears not to 
show NCRI: the two isotopes apparently possess similar 
elastic properties but different inertial properties. This 
isotopic dependence is well documented, in particular by 
the work of West et al. \ This observation would seem 
to also invalidate the present approach. However, the 
tunnelling motion of dislocations is unlikely to proceed 
in a similar manner in the bosonic and fermionic solids. 
In particular, the process of superclimb, ''""'" which may 
assist the formation of connected plastic flow veins, relies 
on the existence of superfluidity in dislocation cores. ■^''*''" 
This mechanism does not operate in solid '^He. 

Kim et al. directly addressed the connexion between 
shear and NCRI in an ingenious experimental arrange- 
ment allowing simultaneous measurements of both quan- 



tities. They observe, in particular, that the response to 
an increase in drive amplitude differs very significantly 
between both properties. However, the drive is not ap- 
plied in an identical manner for both measurements be- 
cause of details of the cell geometry. Soft layers can be 
located at different places and have different conforma- 
tions: they are bound to respond differently. 

Specific experiments can be performed to probe the 
present model. Shear modulus measurements have not 
been performed in a cell geometry for which the plastic 
flow lines would close on themselves in the way they do in 
torsional oscillators. These measurements should reveal 
the existence of supersoft elastic moduli. Equipping a 
torsional oscillator with a floppy membrane as septum to 
interrupt a quantum-coherent flow but not the continuity 
of stresses and strains offers another venue. ' ' The study 
of higher resonance modes in multiply-connected acous- 
tic cavities can also provide a way to probe the internal 
response of an inhomogeneous sample. ' ^ 



Multiple-mode TO resonators' 



appear to give 



somewhat indecisive answers but still show the expected 
trend of enhanced NCRIf at higher frequency. A two- 
mode TO with the dummy massive bob inside the res- 
onator chamber, in contact with the solid helium but 
connected loosely to the main body by an additional tor- 
sion rod provides a mean of coupling shear to the sample 
in a Couette-type experiment. If the inner bob angular 
position could be tracked by some optical or electrody- 
namical means, the coupled system response could be 
analysed in detail. A strong enhancement of G/GcS is 
expected, which would be directly related to the NCRIf. 
A control experiment with bcc '^He, which shows no shear 
modulus anomaly and no NCRI, can be carried out at 
appropriate density and shear modulus values to distin- 
guish between cell and sample contributions to the ap- 
parent NCRI. 

To conclude, the soft layer model presented here takes 
into account known heterogeneities in dislocation pat- 
terns revealed in particular by the anomalous softening 
of most samples of hep ''He. It is argued that the actual 
softness can be even more extreme than observed, being 
hampered by the tortuous arrangement of the dislocation 
structures and of the crystalline regions that they de- 
limit. The corresponding values of the NCRIf are shown 
to lie within the range of observations, barring the high- 
est ones. The model conflicts in no irredeemable way 
with the available assortment of experimental observa- 
tions. Conversely, it can be stated that most existing 
experiments to date support the assumption of the ex- 
istence of very mobile macroscopic veins arranged along 
connected paths in hep crystals of helium 4 and formed 
in a process that depends on quantum statistics, like su- 
perclimb. Matter in the veins themselves undergoes dis- 
placements governed by classical mechanics and subject 
to dissipative mechanisms. These various assertions are 
amenable to experimental veriflcation. 
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Appendix: 

In a crystal lattice with hexagonal symmetry, there 
exists three glide directions for edge dislocations in the 
basal plane, perpendicular to the c-axis, at 27r/3 from 
one another. These dislocations arrange themselves on a 
hexagonal network in the basal plane with side length a, 
possibly connected to adjacent basal planes at a distance 
c along the c-axis by pillars of screw or mixed disloca- 
tions. 

The dislocation network that entirely fills a given basal 
plane of a sample taken as a square of side S for simplic- 
ity, can be constructed as shown in Fig. 5. The building 
block in thicker line is duplicated and translated by A 
along one side and b along the other. There are a total of 
S"^ /bA such translations, each involving a dislocation line 
length equal to 6a. The volume spanned in the process 
is S'^c so that the density of edge dislocation amounts to 
A = 6a/bAc. As b = 3^/^a and A — 3a, there comes that 
Aac = 2/3^/^. This result has already been quoted by 
Iwasa et al. 

The network length Ln can be taken equal to a; it 
disappears in the final result for AL^, which is scale- 
independent. For a hexagonal close packed lattice, c — 
{8/3y^^a so that the expression Ai^ in Eq.(l) takes the 
value 2~^/^. This value is smaller if the lattice is less 
densely packed. For a cubic lattice, a similar derivation 
gives the often quoted geometrical limit AL^ = 3, a value 
much larger than for a hexagonal lattice. This result 
refiects the paucity of easy glide directions in the latter 




ETC. 5. Dislocation network in a hexagonal structure. The 
network is formed of 2D hexagonal cells, of side length a in 
basal planes perpendicular to the c-axis. The whole pattern 
can be generated from the elementary block in thick line of 
length A by translations of moduli A and b as shown in the 
figure. 



